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Abstract 

In this paper we study the problem of calculating the convex hull 
of certain affine algebraic varieties. As we explain, the motivation for 
considering this problem is that certain pure-state measures of quantum 
entanglement, which we call polynomial entanglement measures, can be 
represented as affine algebraic varieties. We consider the evaluation of 
certain mixed-state extensions of these polynomial entanglement mea- 
sures, namely convex and concave roofs. We show that the evaluation 
of a roof-based mixed-state extension is equivalent to calculating a hyper- 
plane which is multiply tangent to the variety in a number of places equal 
to the number of terms in an optimal decomposition for the measure. In 
this way we provide an implicit representation of optimal decompositions 
for mixed-state entanglement measures based on the roof construction. 



1 Introduction 

The burgeoning field of quantum information science is beginning to provide 
us with hints that quantum systems are capable of performing spectacularly 
powerful information processing tasks. An example of this is the existence of 
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a quantum algorithm, Shor's algorithm, which can factor in polynomial time 
Sho94l ISho97| . A consequence of discoveries such as Shor's factoring algorithm 
is the emergence of a widespread belief that, parallel to the physical resource of 
information, there is a corresponding resource which quantum systems support 
called quantum information. The study of how quantum information can be 
manipulated and processed is the central goal of quantum information science. 

So what is quantum information? It is folklore within the quantum infor- 
mation community that it may be quantum entanglement, the unique prop- 
erty possessed by quantum states of bipartite systems not expressible as sim- 
ple tensor products. A classic example of an entangled state is the singlet 
l^ - ) = ^(|01) — 1 10)), which is commonly considered to be the fundamental 
unit of entanglement. There are a number of fascinating information processing 
tasks which use entangled states as a resource, such as teleportation |BBC + 93| 
and enhanced classical communication NCOO], 

Throughout the progress of the research conducted in quantum informa- 
tion science, in particular in the study of quantum entanglement, a number of 
fascinating mathematical problems have been encountered. Within the theory 
of quantum entanglement and quantum communication these problems can be 
broadly described as the the classification of positive operators on tensor prod- 
ucts of vector spaces. This classification problem is central to these theories. 
Many of the stumbling blocks currently encountered by researchers in quantum 
information science are due to incomplete results in this classification program. 

Despite the many problems encountered, a fruitful interplay between math- 
ematics and physics has already resulted from progress in this classification 
program. An example of this is the recent work on indecomposable positive 
linear maps. In this case a result obtained previously by Woronowicz Wor76 
was applied by Pawel Horodecki to construct quantum states which have the 
property that they are bound entangled Hor97 (for a review of bound entangle- 
ment and distillable entanglement see the paper [HH01 ). The subsequent work 
on bound entanglement was then fed back into mathematics by Terhal |Ter01| . 
who used this work in physics to construct new indecomposable positive linear 
maps. 

Bound entanglement is a particular nonlocal phenomena exhibited by certain 
mixed states 1 which falls within the framework of mixed-state entanglement, 
or the study of nonlocal quantum features of mixed states. Understanding 
mixed-state entanglement is an immediate and pressing priority in quantum 
information science. This is because a number of expected developments in 
areas such as quantum communication cannot take place until certain questions 
about mixed-state entanglement are answered. 

Quite apart from the ramifications for quantum information, we believe that 
the study of mixed-state entanglement will lead to new and interesting results 
in mathematics. In particular, there are suggestive indications (sec, for exam- 
ple, the work on concurrence for some striking results |HW971 IWoo98l lUhlOO 




DM01, Osb02j) that there are deep results waiting to be found in the field 

Mixed states are probabilistic mixtures of quantum states. 
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of matrix analysis. 

The purpose of this paper is to provide a method to evaluate certain func- 
tions on the space of positive operators, called mixed-state entanglement mea- 
sures. To do so, in jj^lwe begin by establishing some results in convexity theory. 
These results allow us to state our main result, a necessary condition for con- 
structing convex hulls of certain subsets. In |J3] we make the observation that 
quantum state-space can be represented by the convex hull of a quadratic alge- 
braic variety. Finally, in ^]we combine these two results to provide an implicit 
representation for a family of mixed-state entanglement measures, namely those 
based on the roof construction. 

2 Preliminaries 

In this section we introduce the main object of study for this paper: the convex 
or concave roof. We also establish some general results concerning the construc- 
tion and evaluation of roofs, the most important result being Proposition 12.61 
We illustrate the definitions and results throughout this section in terms of a 
simplified example. 

2.1 Properties of Convex and Concave Roofs 

We begin by establishing some notation. Let V denote the subset of R™ formed 
from the intersection of the zero-sets of a collection of polynomials £q, . . . , £ a £ 
R[xi, . . . , x n ], which we write as 

V = Z(£ ,...,£ a ). (1) 

The set V is known as an affine algebraic variety. We assume that this set is 
compact in the standard topology on R™, irreducible, and nonsingular, i.e., V is 
a nonsingular variety. (For a gentle introduction to varieties and algebraic ge- 
ometry, including definitions of the terms we use here, see |CLQ97| and [Har95 . 
For a more complete reference see [Ha^7|-) We have assumed that V is irre- 
ducible and nonsingular in order to avoid various pathologies associated with 
the definition of the tangent space for these spaces. It should be noted that it is 
possible, with a little extra work, to extend the results of this section to include 
these cases. However, for our applications the assumption of nonsingularity 
imposes no limitations, so we ignore singular varieties. 

Let / be an arbitrary continuous scalar-valued function defined on V. (The 
assumption of continuity is invoked so that subsequent constructions are at least 
continuous.) It suffices for most of our applications to take this function to be a 
polynomial, in which case / is an element of the coordinate ring of V, f £ M.[V]. 
An example variety V and function / is shown in Figure^ The objective of this 
paper is to construct and evaluate a scalar-valued function g which is equal to 
/ on V , but which is defined over all of conv V, the convex hull of V. Obviously 
there are many ways to define such functions. To proceed, we single out two 



3 




X 



Figure 1: The graph (gr/)(r) of f(x, y) — x 3 on Z{x 2 + y 2 — 1). The graph off 
is a space curve in three dimensions. We have indicated this by placing vertical 
lines connecting the set Z[x 2 + y 2 — 1) with the graph. This representation is 
adopted to reinforce the fact that the graph of / is plotted "above" Z(x 2 +y 2 — l). 



special possibilities for g, defined by the formulas: 



and 



(conv/)(r) = Ml ^p j /(x J -) 



(conc/)(r) = sup<^ VV'/fo 



G V, 



i=i 



(2) 



(3) 



where r is an arbitrary point in convU, pj 6 1 is a probability distribution, 
and to is an arbitrary positive integer. A simple consequence of Caratheodory's 
Theorem is that to is bounded above by the dimension of V, to < v = dim V 
|Roc7()| . 



Definition 2.1 ([Uhl98 ). The functions conv / and cone / defined by (0) and 
J2J are called the convex roof and concave roof of /, respectively. We refer to 
the procedure of constructing conv / or cone / as roofing f. 

One of our main objectives is to evaluate conv / or cone / at an arbitrary 
point r £ conv V. Obviously this can be done by simply carrying out the 
optimisation in and (0. Unfortunately, evaluating conv/ at a point r € 
convV according to 10) and 10 is, in general, very difficult. For this reason, in 
this section, we provide an alternative way to evaluate these functions. 
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Figure 2: The three dimensional convex hull conv(gr /) of the graph of fix, y) = 
x 3 on Z(x 2 + y 2 — 1). Note that this is a three-dimensional convex set. The 
upper boundary is shaded red through green and the lower boundary is shaded 
blue. 



We shall spend a lot of time discussing properties of the graph of / on 
denoted gr/, which is the subset of M" +1 defined by 

gr/4{( r ,/( r ))| re V}. (4) 

Associated with gr / is the convex hull of the graph of /, which is the subset of 
M™+ 1 given by 



conv(gr /) = <^ (x, /j) 



3=1 ) 



(5) 



3= 

where Xj € R", Pj is a probability distribution, and we have invoked Caratheodory's 
Theorem to bound the number of terms in the sum by v' + 1 = dim(conv V") + 1. 
See Figure for an example of what conv(grj) looks like for the variety V and 
function / introduced in Figure ^ 

The connection between conv(grj) and conv/ and cone / is provided by the 
following lemma. 

Lemma 2.2. The convex and concave roofs of f are given by 

(conv/)(r) = inf{ju| (r,^) e conv(gr/)} (6) 

and 

(conc/)(r) = sup{/x | (r,/x) G conv(gr/)}, (7) 

respectively. 
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Figure 3: Convex roof for the function f(x,y) = x 3 on Z(x 2 + y 2 - 1). Note 
that the graph of conv/ is the lower boundary of the convex set conv(gr/). 



Proof. The equality of © and J5J) (respectively, (JJJ) and J3J) follows immediately 
from the definition of (0). □ 

Remark 2.3. Lemma 12721 shows us that the convex roof of / is the "lower bound- 
ary" of the convex hull of the graph of /. Similarly, the concave roof cone / is 
the "upper boundary" of conv(gr/). This is illustrated for conv/ in Figure [31 

The convex and concave roofs of / are distinguished from all other functions 
g on conv V equal to / on V by the following result. 

Lemma 2.4 (Lemma A-3, [Uhl98j; [Uhl03 2 ). The convex roof (respec- 
tively, concave roof) of f is, pointwise, the largest convex (respectively, smallest 
concave) function equal to f on V . 

We now turn to the discussion of the main result of this section: an alterna- 
tive way to evaluate conv / and cone /. This result takes the form of a necessary 
condition which limits the space over which the infimum and supremum in J5J 
and © needs to be taken. 

To begin, we make use of the fundamental duality between convex sets and 
intersections of halfspaces to write conv(M), where M C M™ is some bounded 
set, as an intersection in M. n of a collection of halfspaces, 

conv(M)= p| «(„,„), (8) 

2 Armin Uhlmann has proved (private communication) this lemma for all compact V and 
continuous /. 
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where W( m ,w) = {x G W 1 | m • x < lj} is a halfspace, and f2 is the set of all 
halfspaces containing M. Where obvious we abuse notation and write i?( m ,w) G 
£1 to refer to the hyperplane H( m uJ ) = {x G W 1 | m • x = cl>} associated with a 
halfspace W( m , w ) € fi. 

Definition 2.5. A hyperplane 77 in R™ is said to be m-tangent to a nonsingu- 
lar variety (or indeed any sufficiently well-behaved manifold) M C W 1 at points 
X! , . . . , x m 6 M if if contains both the points Xi , . . . , x m G M and their asso- 
ciated tangent spaces T X .M. If m = 2 or 3 then H is said to be bitangent or 
tritangent to M, respectively. 

The set f2 contains three types of halfspace. Firstly, there are those halfs- 
paces which completely contain M and whose associated hyperplane does not 
contain any point of M. Secondly, there are those halfspaces whose associated 
hyperplane contains one or more points of M. These halfspaces are said to 
support M. Finally, of those halfspaces which support M at some collection of 
points xi , . . . , x m there are those whose associated hyperplane is p-tangent to 
M at some subset of p points of {xi, . . . ,x m } (because of the assumption of 
nonsingularity, M has a well-defined tangent space at each point). 

One of our main results, for this subsection, is a necessary condition for the 
hyperplane H associated with some halfspace H G to support M: 

Proposition 2.6. Let H G CI be a hyperplane that supports a nonsingular 
variety (or manifold) M C M" at m points Xi, . . . , x m 6 M. The (translate of 
an) (n — 1)- dimensional vector space defined by H must contain, as subspaces, 
the vector spaces T Xj M, j — l...,m and the vector space formed from the 
vectors Xj — xj., Vj ^ A; = 1, . . . , m. In particular, if H G il supports M at 
Xx, . . . ,x m; it must be m-tangent to M atxi,... ,x m . 

Proof. A hyperplane H G f2 which supports M at m points xi , . . . , x m G M 
contains these points so, in particular, it must contain all the vectors Xj — Xfc, 
7^ k = 1, . . . , to. In addition, if H did not contain T Xj M as a subspace, for 
all j, then some part of H would lie inside the convex hull of M near some point 
Xj, contradicting H G £1. The second statement is simply a restatement of the 
first. □ 

Remark 2.7. Proposition 12.61 provides a necessary condition for a hyperplane 
H G to support M . In general this condition is not sufficient as there exist 
hyperplanes H' S" ft which satisfy the condition of Proposition 12.61 The effect 
of this proposition is to identify a subset of $ C 51 such that the convex hull of 
M can still be found by intersecting members of 

We now specialise the results we have just obtained to describe the structure 
of conv(grj). As demonstrated by Corollary 12. 21 the convex and concave roofs 
of a function / are the top and bottom halves of the boundary of conv(gr/). 
It is useful to obtain a characterisation of this boundary. To do so, we make 
use of the observations above and a lemma of Uhlmann Uhl98 to describe the 
nature of these two hypersurfaces. 



7 



It turns out that the surfaces in M. n+1 defined by either conv/ or cone/ 
are comprised of a union of a family of polytopes (convex hulls of finite sets of 
points) which, pairwise, have zero intersection. To see this 3 , note that any point 
(r,/x) of conv / may be written as a finite sum (r 7/ u) = J2j=l ^i( x ji /( x j))- We 
call this sum an optimal (y + 1)- decomposition of r for conv/. (Where the 
convexity or concavity is obvious we simply refer to this sum as an optimal 
decomposition.) A lemma of Uhlmann (lemma 3, |Uhl98| ) shows us that some- 
thing stronger actually applies, namely that the polytope defined by the convex 
hull of the points (xi, /(xi)), .... (x^+i, f(x v+ i)) also lies entirely within the 
surface defined by conv/. The consequence of this observation is that the 
surface defined by conv/ (which is one half of the boundary of conv(gr/)) is 
made completely from a (pairwise disjoint) union of polytopes. Further, each 
of these polytopes P is a piece of the hyperplane which supports gr / at the ex- 
treme points (xi, /(xi)), . . . , (Xu+i, f(x v+ x)) (because the polytope lives on the 
boundary of the convex set conv(gr/)). The structure of conv / is illustrated 
for the example of / = x 3 in Figure [31 where, in this case, it may be seen that 
conv / is comprised of line segments and one triangle. 

For a given r S conv V , we are interested in determining the polytope P 
on the boundary of conv(gr/) which contains (r, /(r)). Equivalently, given 
r £ conv V we are interested in determining an optimal m-decomposition for 
conv / at r. One way of approaching this problem is to simply search the 
entire space of all sequences of points (xi, /(xi)), . . . , (x„+i, /(x„ + i)) which 
achieve the infimum in (J2J and check to see if their convex hull contains (r, /(r)). 
Proposition ^ . 61 indicates that this is not necessary because we need only search 
through sequences of points (xi,/(xi)), . . . , (x m+ i, /(x m+ i)) for which there 
exists an m-tangent hyperplane passing through them. This statement is the 
content of the following corollary. 

Corollary 2.8. Let r £ conv(T^). // {A^x^} is an optimal m-decomposition 
for conv f at r then there is a hyperplane H passing through Xj , \fj which is 
m-tangent at Xj,Vj. 

2.2 Calculating Optimal Decompositions for Polynomial 
Roofs 

In this subsection we utilise Corollary 12. 81 to provide an implicit representation 
for an optimal decomposition of a point r G convl^ for conv/. We show that 
when / is a polynomial an optimal decomposition for conv / at a point r € 
conv V may be found by simultaneously solving a set of polynomial equations. 
In principle this enables us to invoke the machinery of algebraic geometry, in 
particular Grobner bases methods, to solve the problem CL097, CL098, Har95 
IHa7f7| . 

3 From now on we will phrase most of our results in terms of the convex roof of /. This 
does not entail any loss of generality as all statements about conv / extend in an obvious way 
to cone /. 
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Corollary 12.81 shows that to find an optimal decomposition for conv/ at a 
point r G conv V we need only search over points xi , . . . , x m for which there 
exists an m-tangent hyperplane H passing through these points. If we were able 
to obtain an expression for the basis for the tangent space T x . V at Xj then we 
could reduce this problem to establishing the linear dependence of a collection 
of vectors in the following way. 

Consider the tangent space T X .V at a point Xj € V. Suppose we had a for- 

(i) (?) 

mula, in terms of the components of Xj, for each the basis vectors , . . . , a t 
of T Xj V, where t — dimT Xj V. Using this formula we could rephrase Proposi- 
tion as the statement that for an m-tangent hyperplane H supporting V at 

a) a) 

Xi, . . . ,x m , the matrix R whose rows are formed from the vectors a-f , . . . , a t 
and x 3 — Xfe, Vj 7^ k = 1, . . . , m must have rank less than or equal to n— 1. This 
condition can be imposed by requiring HJ90 that all the (n x n)-submatrices of 
R have zero determinant. If all the entries of R were polynomials 4 in xi, . . . , x m 
then this condition becomes the requirement that a set of polynomials vanish. 

Using the ingredients discussed in the previous paragraph we are now able 
to express an optimal decomposition for a point r € conv V as the simultaneous 
solution to a set of polynomials. The key observation is that it is possible to 
take the components of the basis vectors a%, . . . , at for the tangent space TV 
of an algebraic variety V (embedded in affine space) to be polynomials in the 
coordinate vector x E V. To see this recall |Har95| that the tangent space to 
V , or the Zariski tangent space, is defined to be the null space of the matrix 
Jjk — §^f-, where V = Z(£\, . . . ,£ a ). Because the entries of J are polynomials 
in the components of x, the components of the basis vectors aj for the null space 
of J can be taken to be polynomials in x as well. (To see this, find the basis for 
the null space of J via row reduction. The entries will be rational functions of 
Xj. To get polynomial entries just clear denominators.) 

We now return to the matrix R which we introduced in the previous two 
paragraphs. This matrix can be taken to have m rows consisting of the vectors 5 
Xj — Xj-|_i, j = 1, . . . , m, where we identify j = 1 with j = m + 1, as well as 

mt rows consisting of the basis vectors , . . . , , j — 1, . . . , m. There are 
( m ( i+1 ) ) ( ™ ) n x n-submatrices of R. We denote the determinants of each of 
these submatrices by A a (xi, . . . , x m ), where a runs over all the ( m ^* n f1 ^) (™) 
possible submatrices. In order for a hyperplane to be m-tangent to gr(/) we 
require that all of the polynomials A Q (xi, . . . ,x m ) simultaneously vanish. By 
adjoining the equations tj{~x.k), j = 1, . . . ,a, k = 1, . . . ,m (which ensure that Xj 
lie in V) and the equations = Y^jLi Pj x j — r, = Y^jLi Pj ~ 1 (which ensure 
that r lies in the convex hull of Xj) we obtain a set of conditions for {pj,Xj} to 
be an optimal m-decomposition for r 6 conv V. This discussion constitutes a 
proof of the following proposition. 

Proposition 2.9. An optimal m-decomposition for r S convU consists of m 

4 By which we mean polynomials in the individual components of the vector Xj . 
5 The other vectors Xj — , j '• ^ (k + 1) are linearly dependent with this set. 
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points Xj G 



jn+l 



and m variables pj which satisfy the following equations: 



= ^(x fe ), Vj = 1,. ..,a,Vk = 1, ...,m, 

= A a (x 1 ,...,x m ), Va= 1,..., )(-). 

m 

i=i 

m 

o=£)Pi-i. 



(9) 
(10) 

(11) 
(12) 



Proposition ^. 9l has two important consequences for the evaluation of convex 
roofs. The first is that it provides an implicit representation for the optimal 
decomposition of a point r S convT^ for conv/. This is in contrast to standard 
procedure to calculate an optimal decomposition, i.e., via direct minimisation 
of (J2J. In general this minimisation must be performed numerically and, as a 
consequence, there is no guarantee that a calculated decomposition is optimal 
(It could just be a local minima). The polynomials of Proposition 12 . 91 therefore 
provide a certificate of optimality. 

The second consequence of Proposition 12.91 is that the space required to 
search through to calculate an optimal decomposition is dramatically reduced. 
Often, in practice, an infinite search space is replaced with a finite search space. 

To illustrate Proposition 12 .91 we write out the equations that must be solved 
to find a bitangent plane for the example introduced in Figure ^ These equa- 
tions read 



= ^ + ^-1, 
= ^+^-1, 

yi -xx 

= 2/2 -X 2 

x\ - x 2 yi - y 2 
= pxi + (1 - p)x 2 - 
= pyi + (1 - p)y 2 - ■< 



x\yi 
x\y2 

?1 — 



(13) 
(14) 

(15) 

(16) 
(17) 



By solving these equations we were able to calculate the convex hull and convex 
roof illustrated in Figure [21 and Figure |3J 



3 The Geometry of State Space: the Poincare 
Sphere 

In this paper we are interested in calculating roofs for certain scalar-valued 
functions defined on the state-space of quantum systems. Initially, it seems that 
this problem is unrelated to the results of the previous section because state- 
space (a Hilbert space) has no obvious geometric character. In this section we 
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show that in fact pure state-space for quantum systems of dimension D may 
be represented by the common zero-loci of a collection of polynomial equations. 
This provides the first piece of information — the domain V for the function / 
we wish to roof — needed to apply the results of JSJto our problem. (We discuss 
the second piece of information, the construction of the function / we wish to 
roof, in the next section.) We note that there are many papers dealing with 
generalisations of the Bloch sphere construction for higher-dimensional quantum 
systems. The generalisation we recall here seems to date from 1997, where it 
was considered in detail in the papers of [KMSM97, AMM97]- The observation 
that this embedding of state-space into W 1 is a quadratic variety is, to the best 
of our knowledge, original. 

3.1 Poincare Spheres Based on su(D) 

The Bloch sphere representation for a rank-2 density operator 6 is a fundamental 
conceptual tool in quantum mechanics (see, for example, |NU00l and |Pre 98 for 
a discussion of the Bloch sphere). The fact that the space of all rank-2 density 
operators corresponds to a geometrical object as simple as the sphere S 2 and its 
interior in M 3 has led to many discoveries. This is due, in part, to the fact that 
it is possible to employ geometric intuition when dealing with systems whose 
states or dynamics can be represented on the Bloch sphere. It is therefore very 
desirable to consider generalisations of the Bloch sphere to see if there is a 
simple geometrical object that represents the state-space for a D-level system. 
Unfortunately there is no object which possesses quite the simplicity of the 
Bloch sphere. However, as we argue, there is one natural way to generalise the 
Bloch sphere construction which is useful when discussing certain aspects of 
rank-Z? density operators. 

In this section we discuss a generalisation of the Bloch sphere representation 
for density operators of rank-D, which we call the Poincare sphere (we choose 
this terminology in order to avoid confusion with the Bloch sphere which specif- 
ically refers to S* 2 ). The first step in the Poincare sphere construction is to 
choose an operator basis with respect to which an arbitrary rank-Z) density op- 
erator can be expanded. We choose the set of traceless hermitian generators 7 
A J , j = 1, . . . , D 2 — 1 for the group SU(D). (As we will show, this choice is 
entirely arbitrary. Any other basis for the operator space may be obtained via 
a linear transformation from this one.) We represent the generators in an or- 
thonormal basis {|i> a )}- The generators are given in terms of the set of operators 



6 Any rank-2 density operator p = p|vi)(vi | + (1 — p)\v2) (v2 \ , where \vj) are the eigenvectors 
of p, may be thought of as the state of an effective two-level system (qubit) whose states are 
"logical zero" |0) = \vi) and "logical one" |1) = \v2)- 

7 These matrices are also known as the Gellmann matrices. 




r a ,r; 



+ ) ,r^ ) } defined by 
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ri| } ±-j=(\v a )(v b \ + \v b )(v a \), l<a<b<D, (19) 

^ = -^(\va)(v b \~\v b )(v a \), \<a<b<D. (20) 

The set T consists of D 2 — 1 hermitian operators. We assume that the set T is 
ordered, and we set A 3 to be equal to the jth element of T. Note that, in this 
paper, we suppose that roman indices j run from 1 to D 2 — 1. 

The hermitian operators A J form a representation for the Lie algebra su(D), 
which allows us to write 8 

\3\k = ±. S jk + d jk iX i + ifjk iX l t (21) 

where p k i are the completely antisymmetric structure constants of su(D), and 
the completely symmetric coefficients d? \ are given by 

{A J , A fc } = jj5 jk + 2d jk i\ l . (22) 

By adjoining the operator 

A° = -U/ (23) 



to the D 2 — 1 hermitian generators we obtain a complete basis for the space 
of all D x D hermitian operators. Further, this basis is orthonormal under the 
Hilbcrt-Schmidt inner product: 

tr(A a A /3 ) = S al3 . (24) 

(In the following greek indices are assumed to run from to D 2 — 1.) 

Using the operator basis (|18fl - H20|) we can obtain a geometric representation 
for the space of all pure states of a Z?-dimensional quantum system. Recall 
|Per93| that a hermitian operator p represents a state (either pure or mixed) of 
a quantum system if it satisfies the following conditions: 

tr(p) = l, pi = p, and p 2 < p, (25) 

where by p 2 < p we mean that p — p 2 is positive semidefinite and we require 
equality when p is a pure state. (These conditions assert that p is positive 
semidefinite with trace unity.) Utilising the fact that any density operator may 
be expanded uniquely in terms of the operator basis as 

p =1 Cq A q , (26) 

8 We are employing the Einstein summation convention over repeated indices, so that di k i\ l 
means E^i" 1 d^ M \ l . 
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where the coefficients c„el are given by 

c a = L>tr(pA Q ), (27) 

we can translate the pure-state conditions (|25[1 into geometric conditions on c a . 
Firstly, normalisation implies that 

P=— ^— ■ (28) 

In this way the components c can be regarded as a vector in a (D 2 — 1)- 
dimensional real vector space (this is where we impose the hermiticity require- 
ment), which we can identify with R D _1 . If the density operator is a pure state 
p = \-i/j)(ip\, then tr(p 2 ) = 1 implies that 

\cf=D(D-l). (29) 

Before we complete our description of the geometric conditions we need to 
introduce three vector operations on the vectors c defining pure states. The 
first such operation is simply the euclidean inner product, 

&-h = a] b>. (30) 

The second operation is a generalisation of the cross-product formula in 3 dimen- 
sions, the wedge or outer product, which is constructed using the antisymmetric 
structure constants 

(aAb), = aj ^ ( . (31) 
The third operation is a completely symmetric vector product, the star product 

(a*b)* = aj b k d? k h (32) 

The vector operations (|30|I - H32|I satisfy a number of useful relations. To see this, 
let U £ SU(D) be a group element. Let 0(U) be the operator transforming in 
the adjoint representation of SU{D) that represents U whose matrix elements 
are given by the formula 

[0(U)U = tr(\ a U\W). (33) 

The vector operations defined above transform in the following way under 
the adjoint representation of SU(D), 

0(U)a ■ 0(U)b = a • b, (34) 
0(U)a A 0(U)b = 0(U)(a Ab), (35) 
0{U)a * 0(U)b = 0(U)(a * b). (36) 

By definition the inner product is an SO(D 2 — l)-scalar. Both the star and wedge 
product of two vectors transforming in the adjoint representation of SU (D) give 
rise to another vector transforming in the adjoint representation. 
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We can now state a geometric version of the final pure-state condition. Equa- 
tion 129fl is not a sufficient condition for the vector c to represent a pure state. 
This is a consequence of the fact that the identity tr(p 2 ) = 1 does not imply 
p 2 = p when D > 2. The further constraints due to the anticommutation rela- 
tions (1221) must be taken into account. If these constraints are taken into account 
then the following condition, along with (|29|l . is necessary and sufficient for the 
state represented by c to be pure. The additional constraint on c is known as 
the star-product formula [KMSM971 [AMM97| . 

c*c=(D-2)c. (37) 

With these constraints, the space of all pure states is the subvariety V* C S D _1 
which is the common locus of zeros of the following set of quadratic polynomials, 

£ {c) = |c| 2 - £>(£>- 1), (38) 
£i(c) =c T D l c-efc, (39) 

where D l are real symmetric matrices with entries given by [D ]jk = dp i, and 
ej is a row-vector whose jth entry is one and the rest zero. 

The space V of all pure states may be envisioned as a section V* of the 
sphere S D _1 . The space of all mixed states T>* in this representation (the 
Poincare sphere) is simply given by the convex hull of this quadratic variety, 
V* = convT 5 *. 

To get an idea of what V* looks like as a subspace of S D ~ 2 consider the 
expression 

tr(H = 1 + (40) 

where p — I+ ^' x and a = I+ ^' x are any two pure states. Using the formula 
r ■ s = |r||s| cos 6* we find that the angle 9 between the vectors representing p 
and a is given by 

a (Dtr(pa)-V 



(41) 

D — 1 ' 

This angle is largest when tr(pcr) is zero. In this way we see that widest an- 
gle 6* max between any two points in V* is given by cos -1 (jtijjY so that if r 
represents a pure quantum state, — r cannot, unless D = 2. 



3.2 Poincare Spheres from General Operator Bases 

In the previous subsection we introduced a geometric representation for the 
convex set of all density operators for a D dimensional quantum system. In 
order to do so we expanded an arbitrary density operator p in terms of a specific 
operator basis, the generalised Gellmann matrices. In this subsection we briefly 
summarise the analogous results for when we expand p in terms of an arbitrary 
(hermitian) operator basis. 
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Our main motivation for investigating Poincare sphere representations for p 
with respect to different operator bases is to better study states p of a multipar- 
tite system. Because we are interested in situations when a pure-state density 
operator p is a product state p = a® u>, it is desirable, in some cases, to expand 
p in terms of an operator basis which makes this property manifest, such as the 
basis {X a (£> A' 3 } formed from the tensor product of two Gellmann bases. 

If we expand a density operator p in terms of an operator basis different from 
the Gellmann matrices we find that the generalised Poincare sphere so derived 
is different. It turns out however, that the new Poincare sphere is related to the 
original sphere via an overall isometry. To see this we need to derive a more 
general geometric form of the pure state conditions (|25|l . For clarity we limit 
ourselves to complete hermitian orthonormal operator bases p a , a = 1, . . . , D 2 . 

The geometric pure state conditions analogous to and (|3H follow from 
expanding the (hermitian) product \i a pP in terms of the complete basis, 

MV = X a/ V, (42) 

where are the expansion coefficients. Following the analysis detailed in 

the previous subsection we can write out the pure state conditions in terms of 



c a r a = 1, (43) 
c a c a = 1, (44) 

CaC /3 X Q ' 3 7 = C 7 , (45) 

where t" = tr(^i a ). The pure state conditions (|4~4fl . and (|43|l say that c a 

must be a vector in M. D which lies in the the intersection of the hypersphere 
l|44|) with a hyperplane l|43|l and D 2 quadric surfaces (|45H . 

The algebraic variety V defined by l|43l) . (|44|l . and H45|l representing the 
space of all pure states of a D-dimensional quantum system is related to the 
su(D) Poincare sphere in a simple way. Any complete orthonormal hermitian 
operator basis p a can be found as an orthogonal rotation of the su(D) basis, 

p a = O a p\P, (46) 

where O a p € SO(D 2 ). Say we wish to represent p with respect to some other 
operator basis, i.e. write p — c a p a — d a X a , and we want to know how the 
geometric pure state conditions for c a and d a are related. This relationship can 
be elucidated by making use of the formula c a = d a O a p. By substituting this 
formula into the geometric pure state conditions l|43(l . (|44J) . and Q45JI we find 
that d a must satisfy 

d a r' a = 1, (47) 
d a d a = 1, (48) 

dadpx'^y = d n ( 49 ) 
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where r' Q = O a ^ and x'^ \ = O a a '0^a' a '' 5 'y0 7 ' r (When X a is the 
Gellmann basis these conditions reduce to l|38[) and |39|l.) In this way we see 
that the generalised Poincare representation is unique up to an isomorphism of 
varieties. The geometric representation of the space of all rank- 1? mixed states 
is simply given by the convex hull of the variety V defined by (|T7jl . (j^5|l . and 
JUJJ. We note, in passing, that the variety V defined by is a geometric 
representation of the extreme points of the positive cone of nonnegative D x D 
matrices. The space of all positive linear combinations of points in V is a 
geometric representation for the positive cone of nonnegative D x D matrices. 

4 Entanglement Measures 

In this section we briefly review the definitions and some results that have been 
obtained in the theory of entanglement. In particular we introduce a set of 
axioms for a function to be considered an entanglement measure. We then 
introduce a family of polynomial functions which satisfy these axioms. Finally, 
we demonstrate that the method of S0can be applied to evaluate these functions. 

Suppose we have a bipartite quantum system AB composed of two subsys- 
tems A and B whose respective Hilbert spaces, Ha and He, have dimensions 
dim 7^ = d,A and dim7is = ds- Our problem, at its most elementary level, 
is to answer questions such as: given a pure state p — \ip)(ip\ of AB, when is 
this state expressible as a tensor product p = a (8 cj? More precisely, one of the 
general problems we are interested in is the following. Given an arbitrary mixed 
state p of a multipartite quantum system A1A2 • • • A n , when is p expressible as 
a convex combination of tensor-product states: 

s>rank p 

P= E Pi°a\®°aI®---®°aI (50) 
3=1 

where <t^ are pure states? This problem is fundamental in the theory of entan- 
glement and has received a great deal of attention in recent times. This work 
has culminated in the development of several criteria to determine when a state 
can be written in the form (| 5 f )|> . 

Deciding if a mixed state p is expressible as in ()50|l is significantly compli- 
cated by the fact that there are many ways to write a mixed state as a convex 
combination of pure states. Fortunately there is a recipe which can be utilised 
to relate all the different decompositions of p in terms of convex mixtures of 
pure states |Sch36l |Jay57| II 1.1 \\ !).'<! . To describe this recipe we first write p in 
terms of the decomposition 

r r 

P = J2 Pi\ u i) ( U J I = ^3)^3 1 • ( 51 ) 

where r is the rank of p, {\uj}} are the eigenvectors of p and \uj) = y/pjluj) are 
the subnormalised eigenvectors of p. It turns out that any other decomposition 
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of p can be found via the equation 

P = £to>tol' (52) 
i=i 

where |t7j) = X]fc=i ^jfcfe) an d ^ ^ s an y right-unitary matrix. (The number 
s in Eq. (|52|) may take any integer value greater than or equal to the rank of 
p. This is because p may be realised from an ensemble consisting of arbitrarily 
many pure states.) 

The problem described in the previous paragraphs is known as the separa- 
bility problem — a state p being separable if it is expressible as in (|5t)|l . and 
entangled otherwise. By itself this is an interesting mathematical problem. 
However, there are deeper reasons for why we are interested in answering it. 
These stem from the fact that entanglement, the property possessed by entan- 
gled states, seems to be a fundamentally new kind of resource like free energy or 
information. Motivated by this, researchers are interested in working out when 
a state contains some of this resource (i.e. deciding when p is not separable). 
Ultimately, however, researchers are more interested in working out how much 
of this resource is present in a given state. 

In this area of research, the theory of entanglement, a unified approach to 
answering questions such as the separability of a state p and the measurement of 
the entanglement of p has been developed. This approach is based on studying 
the properties of certain functions, called entanglement measures, on state-space 
which are zero for all separable states, and which are nonzero for all other states. 

At the current time there are only hints that entanglement measures actually 
measure the resource character of entanglement 9 . We ignore this problem for 
the time being, concentrating on the more basic problem of evaluating certain 
entanglement measures. 

4.1 Pure-State Entanglement Measures 

What properties must a scalar-valued function F : Tt — ► M possess in order for 
it to be described as an entanglement measure? This problem has been widely 
studied in the physics literature and a partial consensus has been reached 10 . 

For bipartite quantum systems with state-space isomorphic to C M ®C N , with 
N < M, it has been shown |VidOO| that entanglement measures are in bijective 
correspondence with unitarily invariant, concave (and hence continuous) real- 
valued functions on the space of density operators over C M . More precisely, 
consider any (continuous) real- valued function / : 2?(C M ) — > M on the space 
2?(C M ) of density operators p on C M which satisfies: 

9 This has only been conclusively established for the entanglement for pure-states of a 
bipartite system IBDS W96I IPH97I IvTdOOl INieOOI . 

10 To date, the physics behind the measurement of entanglement for multipartite quantum 
systems has not been satisfactorily elucidated. Only in the bipartite case has agreement been 
reached. The conditions we use for a function to be a multipartite entanglement measure 
follow IVidOOI . 
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1. invariance under the transformation of p by any unitary matrix U : C M — > 
C M , i.e. 

f(UpU^) = f(p), and (53) 

2. concavity, i.e. 

/(p)>p/(wi) + (l-p)/(w 2 ), (54) 
where p = pwi + (1 — p)w2, P € [0, 1], and wi, W2 € 2?(C M ). 

Using / a (pure-state) entanglement measure F on H can be constructed by 
defining 

F(|0» 4/(^(10X01)), (55) 

where tr^ denotes the partial trace over subsystem B. 

Using the correspondence between pure-state bipartite entanglement mea- 
sures and unitarily invariant concave functions we propose a family of entangle- 
ment measures using for / the functions 

/ a (p)^2(l-tr(p a )), (56) 

where a € N, a > 1 . The functions f a are easily seen to be unitarily invariant 
and concave. The family f a has one further property, namely that they are zero 
for pure states p. This means that the entanglement measures F a that they give 
rise to under the construction discussed above are zero if and only if \<f>) is a 
product state. 

If we use the Poincare sphere representation of 7i we see that f a are polyno- 
mial functions of the coordinate vector c a for a state \<j)). This can be established 
as follows. First expand the state \<f>) in terms of the complete orthonormal op- 
erator basis given by {A Q ® X 13 \ a = 1, . . . , M 2 , (3 = 1, . . . , iV 2 }, which is the 
tensor product of the Gellmann matrices of sides M and N respectively: 

p = c ap \ a ® \P. (57) 
Taking the partial trace over B yields 

tr B (p) -c Q0 A Q , (58) 
because A" is traceless for a^O. If we now calculate f a we find that 

/ a (|0))=2(l-tr((c aO A Q r)), (59) 

which is a polynomial in the coefficients c a( g. This function is defined on the 
Poincare sphere for the tensor-product operator basis {A a ® A 73 }. The entangle- 
ment measure f a can also be written as a polynomial in terms of the position 
vector d a defining p in any other operator basis by replacing c a o with Oa^dp, 
where O is the SO{M 2 N 2 ) operator connecting the two operator bases. 
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4.2 Mixed-State Entanglement Measures from the Roof 
Construction 



In the previous subsection we discussed a set of conditions for a function on 
state-space to be considered an entanglement measure. We introduced a family 
of entanglement measures which are polynomial functions of the vector repre- 
senting a state in the Poincare sphere representation. In this subsection we 
discuss a technique to construct entanglement measures defined for all states, 
both pure and mixed, from pure-state entanglement measures. 

Any mixed state p can be written as a convex combination {pj, \ 4>j)} of pure 
states, 

r 

P = I>l&Mil- (60) 

3=1 

As we discussed at the beginning of this section this representation is not unique. 
In order to measure the entanglement of a mixed state p it is necessary to take 
account of this nonuniqucness of the pure state decomposition. Each decom- 
position of p in terms of pure states represents a way of manufacturing the 
mixed state p, because p — Y^ 1 j=iPj\ < t > j)(4'j\ can be thought of as being pro- 
duced from a quantum source emitting the pure state \<pj) with probability 
Pj. If we now imagine that entanglement is an expensive quantity that "costs" 
FQ'fij)) — /(t r s(|0j')(0il))j then it is natural to ask what is the minimal cost of 
producing p. This leads us to define a quantity which is called the (generalised) 
entanglement of formation (see BDSW96 for the original definition and further 
discussions): 

r>s>r 2 

JF(p) =^ inf £ <&/(tTB(|^Wil)), ( 61 ) 
{ijM>j)i ~[ 

where p = X)J=i 9j\' t Pj)(''Pj\j which quantifies the minimum expected cost 11 of 
preparing p. 

If we translate the definition of the mixed state entanglement measure T 
into the Poincare sphere representation we find that it coincides precisely with 
taking a roof of the scalar valued function /. In this way the results of J2Jcan 
be applied to calculate T — conv/. 

5 Conclusions 

In this paper we have considered the problem of calculating the convex hull of 
certain class of affine algebraic varieties. We have shown that this problem is 
intimately related to a problem studied in quantum information science, namely 
that of calculating mixed-state entanglement measures. Our principle result is 
that an optimal decomposition for a point lying in the roof of a function / must 

1 It should be noted that the interpretation of T as quantifying the cost of producing the 
state p is subject to the condition that the entanglement of formation is additive on tensor 
products. 
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satisfy a set of polynomial equations. This provides an implicit representation 
for such decompositions. 

It is tempting to apply the well-developed techniques in computational com- 
mutative algebra to calculate the polynomial entanglement measures we intro- 
duced in 21 for mixed states. Unfortunately, apart from the example we consid- 
ered in Figures n 13 and |3 (which pertains to a cubic entanglement measure for 
real rank-2 density operators) the polynomials that appear in Proposition 12.91 
are far too complex for current computational approaches. 

One principle future direction suggests itself at this stage. Determining 
the additivity of convex functionals (particularly entanglement of formation-like 
quantities) on tensor products of quantum states is a major unsolved problem 
in quantum information science. It is plausible that the implicit representation 
afforded by Proposition 12 . 91 may provide one route to considering this problem. 
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